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Effects  of  a Time -Varying  Photoelectron 
Flux  on  Spacecraft  Potential 


1.  INTRODUCTION 

Although  the  potential  between  a satellite  and  the  ambient  plasma  is  usually  on 
the  order  of  a few  volts,  the  potential  can  increase  to  values  as  high  as  -18,000  V 
on  time  scales  of  a few  minutes  as  the  satellite  passes  in  and  out  of  the  earth's 
shadow.  During  such  periods,  increased  arcing  occurs.  * For  this  reason,  an 
understanding  of  the  processes  involved  in  these  rapid  potential  variations  is  of  real 
and  immediate  interest  to  spacecraft  designers.  The  purpose  of  this  report  is  to 
present  a model  of  these  variations. 

In  order  to  model  the  potential  variations  experienced  by  a spacecraft  as  it 
enters  or  exits  the  earth's  shadow,  it  is  necessary  to  estimate  the  changes  in  the 
solar  illumination  of  the  spacecraft.  This  allows  the  determination  of  the  photo- 
electron  current  variations  and  the  resulting  potential  variations.  In  the  first  part 
of  the  report  a model  of  the  atmospheric  attenuation  of  the  solar  illumination  is 
developed.  In  a previous  study  a model  capable  of  estimating  satellite  potentials 
given  the  ambient  and  photoelectron  currents  was  derived.  This  model  is  combined 
with  the  attenuation  model  to  estimate  the  varying  spacecraft  potential.  The  pre- 
dicted potentials  are  then  compared  with  the  actual  data.  The  results  of  this 

(Received  for  publication  10  May  1978) 

1.  Shaw,  R.R.,  Nanevicz,  S.  E. , and  Adamo,  R.C.  (197  6)  Observations  of  elec- 
trical discharges  caused  by  differential  satellite  charging,  AIAA  Progress  in 
Astronautics  and  Aeronautics,  A.  Rosen  (ed.  ),  47:61-7  6. 


comparison  indicate  a standard  deviation  of  approximately  ± 700  V between  the 
predicted  and  observed  values.  The  FORTRAN  programs  for  the  model  are  included 

In  Appendix  A. 


2.  SATELLITE  OBSERVATIONS 

Three  different  satellites  served  as  data  sources  for  the  study.  The  geosyn- 
chronous satellites  ATS-5  and  ATS-6  were  the  principal  sources  for  data  on  the 
ambient  plasma  environment  and  on  the  satellite  potential,  while  the  polar-orbiting, 

low-altitude  INJUN  5 satellite  provided  information  on  the  photoelectron  current. 

2 1 
The  reader  is  referred  to  DeForest  and  Mcllwain  for  details  of  the  ATS-5  instru-  < 

3 4 

ment,  Mauk  and  Mcllwain  for  the  ATS-6  instrument,  and  Burke  et  al  for  the 

INJUN  V instrument. 

The  ATS-5  and  ATS-6  satellites  are  in  low  inclination,  geosynchronous  orbits. 

( 

The  ATS-5  satellite  is  cylindrical  (1.8  m long  by  1.5  m in  diameter).  It  spins  at  1 

1.  27  rps  with  its  axis  parallel  to  the  earth's  rotation  axis.  ATS-6,  in  contrast,  is 

essentially  a large  (10  m)  spin-stabilized  dish  antenna.  The  plasma  data  were  ob- 
tained from  the  University  of  California  at  San  Diego  (UCSD)  plasma  experiments  ; 

on  ATS-5  and  ATS-6.  Both  instruments  consist  of  electrostatic  analyzers  designed 
to  measure  the  positive  ion  and  electron  populations  between  51  eV  and  51  keV  for 
ATS-5  and  between  0 and  80  keV  for  ATS-6. 

The  ATS-5  and  ATS-6  instruments  measure  the  particle  flux  in  64  channels 
(2  background  channels  and  62  energy  channels),  returning  a complete  energy 
spectrum  in  roughly  15  to  20  seconds.  For  ATS-5,  the  center  energy  for  a channel 
is  112  percent  of  the  center  energy  of  the  previous  channel.  For  ATS-6  this  value 
is  ~ 1 1 3 percent.  This  results  in  an  uncertainty  of  about  ±5  percent  in  the  energy 
(or  potential)  determination. 

The  ATS-5  and  ATS-6  instruments  return  count  rate  data  which  are  converted 
to  differential  energy  spectra.  From  the  differential  energy  spectra,  the  ambient 
currents  and  distribution  functions  necessary  in  estimating  the  satellite  potential  are 
determined.  The  satellite  spectra  also  indicate  the  spacecraft  potential.  Briefly, 
the  low  energy  ion  (electron)  population  is  accelerated  if  the  satellite  has  a negative 
(positive)  potential  relative  to  the  ambient  plasma.  This  acceleration  produces  a 
pronounced  peak  in  the  low  energy  ion  (electron)  channels  at  an  energy  in  electron 
volts  (eV)  corresponding  to  the  satellite  potential  in  volts  (V).  Thus,  ATS-5  and 

2.  DeForest,  S.  E.  , and  Mcllwain,  C.  E.  (1971)  Plasma  clouds  in  the  magnetosphere, 

J.  Geophys.  Res.,  76(No.  16):3587  -361 1. 

3.  Mauk,  B.  H.  , and  Mcllwain,  C.  E.  (1975)  ATS-6  UCSD  Auroral  Particles  Experi- 

ment, IEEE  Trans.  Aerospace  and  Electronic  Systems,  AEA- 1 l(No.  6); 

ii25-n7Tn  

4.  Burke,  W.  J. , Donatelli,  D.  E. , and  Sagalyn,  R.  C.  (1978)  INJUN  5 observations 

of  low-energy  plasma  in  the  high-latitude  topside  ionosphere,  to  appear  in 
J.  Geophys.  Res. 


8 


ATS -6  provide  information  on  both  the  ambient  plasma  and  the  spacecraft  potential. 
The  25  eclipse  passages  for  which  these  quantities  were  determined  are  listed  in 
Table  1. 


Table  1.  Date  and  Approximate  Start  and  End  Time  (UT) 
of  ATS -5  and  ATS -6  Eclipses.  Potentials  immediately 
following  entry  and  preceding  exit  from  eclipse 


ATS-5 

Eclipse  Passages 

Date 

Entry 

Potential 

Exit 

Potential 

22  Sep  1969 

0628 

-3404 

0732 

-3810 

16  Oct  1969 

0626 

-5357 

07  12 

-3810 

12  Sep  1970 

0630 

-2423 

07  19 

-1730 

15  Sep  1970 

0625 

-877 

0722 

-1539 

17  Sep  197  0 

0622 

-537  8 

0723 

-3041 

19  Sep  1970 

0620 

-2725 

07  24 

-1939 

26  Sep  1970 

0725 

-2168 

5 Oct  1970 

0615 

-5357 

0719 

-3811 

17  Oct  1970 

0628 

-1232 

0700 

-2172 

18  Oct  1970 

0632 

-397 

0655 

-316 

19  Oct  197  0 

0638 

-558 

0649 

-558 

ATS-6 

i 

28  Feb  197  6 

2143 

-9330 

2200 

-8140 

6 Mar  1976 

2124 

-7  110 

2219 

-9330 

The  INJUN  5 satellite  is  a magnetically -aligned,  polar -orbiting  satellite  with 
a spin  period  of  20  min,  an  inclination  of  81®,  an  apogee  of  2453  km,  and  a perigee 
of  677  km.  Data  were  used  from  the  Air  Force  Geophysics  Laboratory  (AFGL) 
spherical  electrostatic  analyzers  (SEA).  These  analyzers  consist  of  a spherical 
tungsten  collector  surrounded  by  two  concentric  wire  mesh  grids.  The  current  to 
the  collector  grid  is  measured  at  3 or  40  samples  per  second.  When  the  ion 
collector  is  illuminated,  the  resulting  photoelectron  flux  is  the  dominant  source  of 
the  measured  current.  As  the  spacecraft  approaches  the  earth's  shadow,  the 
photoelectron  flux  decreases  in  response  to  atmospheric  attenuation.  As  will  be 
discussed,  this  will  allow  a calibration  of  the  atmospheric  attenuation  model. 
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3.  ATMOSPHERIC  ATTENUATION 


Theory  and  observation  indicate  a linear  relation  between  the  photon  flux  falling 
on  a surface  and  the  photoelectron  flux  emitted  by  that  surface.  As  a satellite  passes 
into  eclipse,  the  percentage  of  the  solar  disk  that  is  visible  decreases.  The  per- 
centage decrease  in  the  visible  solar  disk  is  equal  to  the  percentage  change  in  the 
photon  flux.  Given  an  estimate  of  the  percent  attenuation  of  the  solar  disk  and  of 
the  total  photoelectron  current  for  the  unattenuated  solar  disk,  it  is  straightforward 
to  calculate  the  changing  photoelectron  currents.  In  this  section  we  estimate  the 
attenuation  of  the  solar  disk.  In  Section  5 we  will  determine  the  total  photoelectron 
current. 

The  problem  of  finding  the  percentage  of  the  solar  disk  that  is  obscured  as  the 
satellite  passes  into  the  earth's  shadow  has  two  parts— determination  of  the  position 
of  the  center  of  the  solar  disk  and  determination  of  the  atmospheric  attenuation.  The 
two  operations  can  be  treated  separately  if  the  concept  of  X , the  minimum  ray 
path  altitude,  is  introduced.  (Figure  1)  is  defined  as  the  minimum  distance  from 

the  earth  that  a ray  of  light  travels  in  going  from  the  center  of  the  sun  to  the  satel- 
lite. The  problem  of  finding  the  percentage  of  the  solar  disk  that  is  obscured  re- 
duces to  that  of  finding  the  atmospheric  attenuation  as  a function  of  Xm  and  finding 
Xm>  as  seen  from  the  satellite,  as  a function  of  time  or  position  in  orbit. 
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Figure  1.  Illustration  of  the  Meaning  of  x»  the  Angle  Between  the  Sun  (S)  and  the 
Earth's  (E)  Center,  Xs,  the  Distance  From  the  Satellite  (SAT)  to  the  Center  of 

the  Earth,  and  X , the  Minimum  Ray  Path 
m J 


l 


T 
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3.1  Determination  of  Xn| 

X is  determined  as  a function  of  satellite-earth-sun  position  in  spherical 
m 

coordinates.  Step  one  is  the  determination  of  the  angle,  X,  between  the  center  of 
the  earth  and  the  sun-satellite  line  (see  Figure  2).  x is  a function  of  the  latitude  of 
the  sun  (6^)  and  the  satellite  (6^1  relative  to  the  earth  and  the  longitude  of  the  sun 
(«q)  and  the  satellite  (oO  (Figure  3).  The  sun  is  considered  to  be  at  infinity  so  that 
the  longitude  and  latitude  of  the  sun  relative  to  a coordinate  system  based  on  the 
earth-mean  sun  (or  satellite-mean  sun)  line  and  the  spin  axis  of  the  earth  are  equal 
for  the  earth  and  satellite.  Thus,  the  latitude  of  the  sun  in  the  satellite  coordinate 
system  is  (Figure  3a) 

S^sat)  = 6^  (earth)  . (1) 


Figure  2.  Detail  of  the  Measurement  of  the  Angle  x 
(S  - SAT  - E) 


The  latitude  of  the  earth  (6^)  in  the  satellite  system  becomes  (Figure  3b) 

6@(sat)  * -6s  (earth)  , (2) 
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Figure  3a.  Meridian  Projection  of  the  Satellite-Sun- 
Earth  System  Showing  the  Measurement  of  , the 
Solar  Latitude,  and  6 , the  Satellite  Latitude 


Figure  3b.  Polar  Projection 
of  the  Satellite-Sun-Earth 
System  Showing  the  Measure  - 
ment  of  o/q  , the  Longitude 
of  the  Sun,  and  a , tne 
Longitude  of  the  Satellite 
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I 

1 
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It  longitude  is  reckoned  in  the  name  direction  from  mean  local  noon  in  the  satellite 
and  earth  systems  (remember,  as  the  sun  is  at  infinity,  the  satellite-mean  sun 
and  earth-mean  satellite  sun  lines  are  parallel),  then  the  longitude  of  the  sun  in  the 
satellite  system  is 

«Q(sat)  * Oq  (earth)  . (3) 

Likewise,  the  longitude  of  the  earth's  renter  in  the  satellite  coordinate  system  is 

o'  (sat)  n (earth)  -180"  . (4) 

(D  s 

Thus,  the  separation  in  longitude  of  the  sun's  renter  and  the  earth's  renter  as  seen 
from  the  satellite  is 


n (Y_(sat)  - om(sat) 

(J  O 

(nt  (earth)  -180°)  - ou  (earth)  . (3) 

s fc) 

We  now  have  the  latitudes  of  the  sun  and  earth  in  terms  of  the  satellite  coordinate 
system  and  their  relative  displacement  in  longitude.  Referring  to  Figure  4,  the 
problem  reduces  to  solving  a spherical  triangle,  the  solution  of  which  is 


cos  x cos  (90’  - 6^)  cos  (90"  - 6^) 

+ sin  (90"  - 60)  sin  (90’  - cos  n 

sin  (6C0)  sin  (6ft)  + cos  (6Q)  cos  (6^)  cos  n . 

Once  x is  given,  follows  immediately  (Figure  1) 

X - XS  • sin  (X>  -R~  , 
m fc, 

where 

XS  «■  distance  from  earth  to  satellite  (in  km)  , 

Hj,  * radius  of  earth  (6378  km)  . 


(6) 


(7) 


r 


13 


Figure  4.  Illustration  of  the  Actual 
Spherical  Trigonometric  Problem  in 
Terms  of  the  Satellite  Reference  System 


t 


1 


X2  Determination  of  Atmospheric  Attenuation 

The  determination  of  atmospheric  attenuation  is  uncertain  since  ihe  exact  re- 
sponse of  satellite  materials  to  solar  flux  is  not  well  known  so  that  the  wavelength 
ranges  of  interest  are  not  well  defined.  Also,  the  apparent  size  of  the  sun  varies 
slightly  with  wavelength  and  with  solar  activity  (prominences  and  the  corona  can 
introduce  asymmetries  in  the  shape  of  the  sun  at  ultraviolet  wavelengths).  Further, 
the  atmosphere  itself  varies  in  time.  The  best  solution  at  present  is  an 
iterative  process  in  which  each  of  the  unknowns  is  varied  within  a reasonable 
range. 

The  process  of  determining  atmospheric  attenuation  has  three  steps.  First, 
a model  of  the  attenuation  of  the  atmosphere  is  assumed.  Secondly,  the  attenuation 
factor  is  integrated  over  the  disk  of  the  sun.  This  is  done  for  the  required  values 
of  X to  give  a table  of  percent  illumination  (assumed  to  be  directly  proportional 
to  photoelectron  current)  vs  X^,  In  the  third  step,  to  be  carried  out  in  a later 
section,  the  results  are  compared  with  a plot  of  the  photoelectron  current  required 
to  give  the  observed  potentials  as  a function  of  X . The  atmospheric  model  and 
total  photoelectron  current  are  varied  until  agreement  is  reached. 

The  satellites'  response  curves  are  assumed  to  be  similar  to  that  of  aluminum 
so  that  only  the  wavelength  range  between  ~ 1000  i and  -.1000  i is  significant.  This 
means  only  ultraviolet  radiation  will  be  considered.  The  cross-section  for  absorp- 
tion is,  in  general,  orders  of  magnitude  larger  than  scattering  at  these  wavelengths 
so  the  radiation  is  attenuated  long  before  significant  scattering  or  refraction  takes 
place.  Likewise,  in  the  altitude  range  over  which  significant  attenuation  is  ex- 
pected, 0^  and  O.j  are  the  dominant  absorbing  species.  The  altitude  (30  to  90  km) 
where  these  species  absorb  most  of  the  flux  is  little  affected  by  latitude  or  geo- 
magnetic effects. 
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I or  an  atmosphere  which  falls  off  exponentially  in  altitude,  the  photon  flux  0 
as  a function  of  wavelength  and  X for  a single  < onstituent  is 

0 (A,  X^l  - 0^  < A)  e"  7 (« ) 


I 


< 


where 


T = n(X  ) a (277  R H) 
m a e 


1/2 


R =■  earth's  radius  , 

e 

n (X  ) - number  density, 
m •' 

II  - scaie  height  of  constituent  at  X , 

m 

0^1  A)  unattenuated  solar  flux  at  wavelength  A , 


For  simplicity,  an  exponential  atmosphere  is  assumed 


n (X  ) = n e 
m o 


-(X  - X ) 

m o 

TX — 


n , X , 6 X constants  for  given  constituent, 
o o n 


Then 

-(X  - Z ) 

-e  m ° 

0 (X^,  A)  - 0w  (A)  e ^ (9) 

where  all  the  constants  have  been  lumped  into  Z and  6 Z. 

Veither  the  0 2 or  0^  density  follows  a strictly  exponential  curve  in  the  region 
of  interest.  The  apparent  diameter  of  the  sun,  however,  is  370  km  as  viewed  by 
a satellite  at  geosynchronous  orbit.  Thus,  a gross  averaging  will  take  olace  such 
that  wavelength  and  scale  height  variations  are  smoothed  out  if  the  attenuation  is 
considered  for  the  entire  solar  disk.  To  first  order,  a model  of  atmospheric  atten- 
uation of  the  functional  form  given  in  Eq.  (9)  should  be  adequate. 

Once  the  attenuation  function  in  terms  of  > is  assumed,  an  integration  over 
the  disk  of  the  sun  can  be  performed.  Again,  .he  calculation  is  carried  out  in 
spherical  coordinates.  Previously  the  angular  separation  between  the  center 

of  the  solar  disk  (assumed  to  be  the  optical  center)  and  the  center  of  the  earth,  was 

5.  Swider,  W.  (1964)  The  determination  of  the  optical  depth  at  large  solar  zenith 
distance.  Planet.  Space  Sci.  12:761-782. 

■ SAAto 
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X = XS  • sin  (X1)  -Re  . 

XI  = cos  _1  (cos  (9)  cob  (x)  + sin  ( 9 ) sin  (x>  cos  («3>)  , 

« ^ = angular  radius  of  sun  (function  of  day  of  year)  , 
r w 

0^,  Zq,  67,  assumed  constants  for  atmospheric  model  . 


The  geometry  is  illustrated  in  Figure  5.  The  percentage  of  solar  flux,  P (X  ), 
given  by  dividing  (X  ) by  0T  (<«).  As  <J>'  (9  ,<p)  goes  to  <PM  as  Xm  goes  to  infinity, 
this  gives 


p«m’  ■ ffwr 


>T  IXn,’ 
l it  1 1 -cos!  a 


(11) 


Figure  5.  Geometric  Representation  of  the  Obscurations  of  the  Solar 
Disk  Upon  Eclipse  Entry.  The  variables  are  defined  in  the  text 
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As  discussed  earlier,  INJUN  V data  provide  an  estimate  of  PIX^I.  In  Fig- 
ure 6,  P (X  ) is  plotted  for  the  INJUN  V satellite  for  the  case  of  no  atmosphere 

(complete  attenuation  below  X =0  km),  for  an  atmospheric  model  having 

m c 

Zq  = 90  km  and  6 Z = 40  km,  and  for  an  0^  atmosphere  (6  Z -4,67  km,  Zq  -7  6.4), 

Also  shown  is  the  percentage  of  photoelectron  current  measured  as  the  satellite  was 

eclipsed.  Even  though  INJUN  5 is  at  2500  km  altitude  where  the  attenuation  should 

readily  show  the  effects  of  atmospheric  variations,  the  simple  model  (assuming 


% ILLUMINATION 

0 20  40  60  80  100 


Figure  6.  Percentage  of  the  Solar  Illumination  Reaching 
the  Spacecraft  and  the  Percentage  of  Photoelectron  Current 
Measured  by  INJUN  5 as  Functions  of  X , the  Minimum 
Ray  Path  Altitude 


6.  Weeks,  L.  H. , Cuikay,  R.S,,  and  Corbin,  J.  R.  (1972)  Ozone  measurements  in 

mesosphere  during  the  solar  proton  event  of  2 November  1969,  J.  Atmos.  Sci.  , 
j!!|<No.  6):  1 138-1 142. 
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7.  - 90  km  and  5 Z = 40  km)  is  adequate,  demonstrating  the  smoothing  over  wave- 

length and  height.  In  lieu  of  knowledge  of  the  photo -emission  response  functions 
for  the  ATS-5  and  ATS-6  satellites,  the  INJUN  5 (or  tungsten)  attenuation  function 
(6  Z = 40  km  and  ZQ  = 90  km)  will  be  assumed.  We  will,  however,  compare  the  re- 
sults for  this  assumption  with  those  for  no  atmosphere.  At  geosynchronous  orbit, 
the  difference  between  the  two  assumptions  is  not  as  great  as  implied  by  Figure  6. 
Also,  the  INJUN  5 attenuation  function  gives  adequate  results  over  the  X range  of 
interest  and,  unlike  other  functions  that  might  be  assumed,  is  based  on  observational 
data. 


4.  SPACECRAFT  CHARGING  MODEL 


Given  the  atmospheric  attenuation  function  we  require  a model  that  predicts 
potential  variations  in  the  presence  of  photoelectron  currents.  This  model  will  be 
needed  in  calculating  the  residual  current— presumably  the  photoelectron  current— as 
the  spacecraft  passes  in  and  out  of  eclipse.  The  total  photoelectron  current  will  be 
determined  by  fitting  the  attenuation  function  to  these  residuals.  Given  this  value 
and  the  attenuation  function,  we  may  estimate  the  potential  variations  as  the  space- 
craft passes  in  and  out  of  eclipse— the  goal  of  the  study.  A model  capable  of  calcu- 
lating spacecraft  potentials  as  a function  of  the  ambient  currents  and  meeting  these 

7 

objectives  is  briefly  presented  in  this  section.  The  reader  is  referred  to  Garrett 
for  more  details. 

The  problem  <■»*  determining  the  potential  between  a spacecraft  and  the  ambient 
plasma  is  conceptually  simple.  The  basic  problem  is  to  find,  for  a given  set  of 
ambient  plasma  conditions,  a satellite  potential  for  which  there  is  0 net  current  to 
the  spacecraft  surface.  Specifically,  it  is  necessary  to  find  V such  that 


J (V)  - (J.(V)  + .1  (V)  + J^.m  * Jot;<V)  + J . (V))  = 0 , 

e I se  SI  BS  ph 


where 


(12) 


J incident  electron  current, 

e 

Jj  = incident  ion  current, 

Jge  - secondary  electron  current  emitted  by  incident  electrons, 

Jgj  " secondary  electron  current  emitted  by  incident  ions, 

backscattered  electron  current, 

* photoelectron  current. 


J 


BS 


J 


ph 


7.  Garrett,  H.  B.  (1978)  Spacecraft  Potential  Calculations  - A Model,  AFGL-TR- 
78-0116.  
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The  currents  are  all  functions  of  the  ambient  particle  distributions  and  satellite 

potential,  V,  Their  calculation,  in  terms  of  the  ATS-5  and  ATS-6  UCSD  plasma 

7 

data,  is  described  in  Garrett.  Given  these  currents,  the  potential  is  determined 
by  an  iterative  procedure  in  which  V is  varied  until  Eq.  (12)  holds. 

The  ratios  of  the  secondary  and  backscattered  currents  to  the  ambient  current 
are  not  known  accurately.  Likewise,  the  photoelectron  current  is  not  known.  It  is 
necessary,  therefore,  to  initially  "calibrate"  the  potential  model.  This  is  done  by 
first  choosing  particle  spectra  immediately  before  (immediately  after)  and  immedi- 
ately after  (immediately  before)  the  spacecraft  pass  into  (out  of)  the  earth's  shadow 
(see  Table  1).  The  photoelectron  current  is  assumed  to  be  0 and  the  secondary 
and  backscattered  ratios  adjusted  until  the  predicted  potentials  agree  with  the  ob- 
served potentials  in  a least  squares  sense.  The  results  of  this  procedure  for  the 
ATS-5  and  ATS-6  satellites  are  presented  in  Figure  7. 


Figure  7.  The  Predicted  and 
Observed  Potentials  for  ATS-5 
and  ATS-6  for  Total  Eclipse  of 
the  Sun.  The  closed  symbols 
are  for  values  calculated  using 
the  eclipsed  spectra.  Open 
symbols  are  for  values  calculated 
using  the  sunlit  spectra 


The  error  bars  in  Figure  7 correspond  to  a ± lo  deviation.  For  the  entire  data 
set  (both  sunlit  and  eclipsed  spectra),  the  error  bar  is  t 1500  volts.  Ignoring  the 
two  outlying  points  (the  plasma  changed  significantly  during  eclipse  passage  for  these 
points),  the  error  is  ± 1300  V for  the  sunlit  spectra  and  ± 5 00  \ for  the  eclipsed 
spectra. 

The  difference  between  the  predictions  due  to  the  sunlit  and  eclipsed  spectra  is 
believed  to  result  from  the  lack  of  knowledge  about  the  particles  outside  the  detector's 
energy  range  and  from  the  finite  interval  (particularly  at  high  energies!  between 
energy  steps.  Shifting  the  sunlit  spectra  in  energy  to  account  for  the  satellite  poten- 
tial accentuates  these  problems.  A similar  difficulty  arises  when  the  eclipsed 


19 


spectra  are  shifted  down  in  potential.  Thus,  an  error  of  about  25  percent  is  intro- 
duced that  is  dependent  on  satellite  potential.  Although  spectra  taken  at  the  same 
time  as  the  desired  potential  calculation  are  clearly  to  be  preferred,  it  is  not 
possible  to  do  this  in  most  cases.  In  the  discussions  that  follow,  the  results  pre- 
sented will  either  be  the  average  of  the  eclipsed  and  sunlit  spectra  or  designated  as 
sunlit  and  eclipsed  estimates. 

In  summary,  given  the  ambient  plasma  distribution,  the  currents  to  the  space- 
craft can  be  estimated  as  a function  of  satellite  potential.  Given  the  currents  as 
functions  of  spacecraft  potential,  the  spacecraft  potential  can  be  calculated.  It  is 
necessary,  however,  to  "calibrate”  the  secondary  and  backscattered  currents.  This 
was  done  for  eclipse  conditions  when  the  photoelectron  current  was  zero.  The  model 
can  then  be  used  to  estimate  the  residual  (photoelectron)  currents  as  the  spacecraft 
is  eclipsed  or  to  estimate  satellite  potentials. 


5.  RESULTS 

As  discussed  in  the  previous  sections,  we  need  to  determine  the  total  ATS- 5 

and  ATS- 6 photoelectron  currents.  Once  given  these  currents,  the  atmospheric 

attenuation  function  gives  the  photoelectron  current  as  a function  of  X . Knowing 

the  ambient  conditions,  we  then  calculate  the  spacecraft  potential  as  a function  of 

X . 
m 

To  determine  the  total  ATS-5  and  ATS-6  photoelectron  currents,  the  spacecraft 
charging  model  is  used  to  estimate  the  residual  current  — that  is,  the  current  neces- 
sary to  give  the  observed  potential  after  all  known  currents  are  accounted  for  — as  a 
function  of  X . In  Figures  8a  and  8b  we  have  plotted  these  currents  (actually  the 
number  fluxes).  Shown  for  both  satellites  are  the  INJUN  V atmospheric  attenuation 
function  and  the  attenuation  function  for  no  atmosphere.  The  attenuation  functions 
have  been  fit  to  the  data.  In  the  case  of  ATS-5,  the  INJUN  5 function  fits  the  ob- 
served fluxes  quite  well,  while  for  no  atmosphere  the  function  does  not  fit  values 
having  an  of  0 or  greater.  For  ATS-6,  it  is  not  possible  to  unambiguously  choose 

either  function.  The  total  photoelectron  fluxes  implied  by  these  fits  are,  in  number/ 

2 , 
cm  -sec-ster 


No  Atmosphere 
(Zq  = 90  km, 


ATS-5  ATS-6 

5X  10‘  1.  2 X 108 

2.  5X  108  4.  OX  108 


Z - 40  km) 


Figure  8a.  The  Residual 
(Photoelectron)  Current 
Estimated  for  21  ATS-5  Eclipse 
Passages.  This  is  the  current 
necessary  to  explain  the 
observed  satellite  potential 
after  all  ambient  and  secondary 
currents  have  been  subtracted. 
Two  atmospheric  attenuation 
profiles  are  also  plotted 


Figure  8b.  The  Residual 
(Photoelectron)  Current 
Estimated  for  4 ATS-6  Eclipse 
Passages.  This  is  the  current 
necessary  to  explain  the 
observed  satellite  potential 
after  all  ambient  and  secondary 
currents  have  been  subtracted. 
Two  atmospheric  attenuation 
profiles  are  also  plotted 


For  ATS-5.  these  values  can  be  converted  directly  to  currents  by  multiplying 

5 2 

by  the  total  spacecraft  area  (1,8  / 10'  cm  ).  dividing  by  the  sunlit  cross-section 
4 2 

(~2.7  / 10  cm  for  the  side),  and  multiplying  by  ffq  (q  is  charge).  The  ATS-5 

2 2 8 

fluxes  give  0.  552  nA/cm  and  0.  11  nA/cm  . Grard  et  al  predict  currents  on  the 
2 2 

order  of  0.4  nA/cm  (graphite!  to  4.  2 nA/cm  (aluminum  oxide).  Thus  the  INJUN  5 
attenuation  factor  gives  low  but  reasonable  values  for  the  total  photoelectron  current. 
Considering  the  complex  ATS-6  geometry,  it  is  not  possible  to  accurately  calculate 
the  current. 

8.  Grard,  H.J.L.  , Knott,  K.  , and  Pedersen,  A.  (1973)  The  influence  of  photo- 
electron and  secondary  electron  emission  on  electric  field  measurements 
in  the  magnetosphere  and  solar  wind.  Photon  and  Particle  Interactions  with 
Surfaces  in  Space,  H.J.L.  Grard  (ed. ),  pp  163-189. 
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The  major  source  of  errors  in  the  preceding  calculation  is  the  satellite  orbit. 
Although  the  calculation  of  the  orbit  of  a geosynchronous  satellite  is  straightforward 
since  relatively  few  perturbations  occur,  the  satellite  ephemerides  are  only  updated 
on  a 2 or  3 week  basis.  All  of  the  orbital  periods  studied  were  4-  to  7 -day  extrapo- 
lations. Errors,  particularly  in  the  time  of  local  midnight,  on  the  order  of  ~ 10  sec 
were  typical  for  ATS -5,  while  the  error  in  the  ATS -6  ephemerides  were  on  the 
order  of  90  seconds.  Such  errors  translate  to  errors  on  the  order  of  ± 15  km  and 
± 150  km  in  Xm>  It  was  necessary  to  assume  eclipse  symmetry  around  midnight 
and  adjust  the  time  of  local  midnight  accordingly.  Although  this  significantly  im- 
proved individual  eclipse  pairs  (particularly  for  ATS-6),  errors  in  latitude  of  ± 0.  05° 
(corresponding  to  errors  in  Xm  of  35  km  or  less)  remain. 

Given  the  INJUN  5 attenuation  function  and  the  corresponding  total  photoelectron 
current  we  use  Eq.  (12)  to  calculate  the  spacecraft  potential  as  a function  of  the 
ambient  plasma  conditions. 

Figures  9,  10,  and  11  represent  the  end  product  of  this  analysis.  In  Figure  9, 
the  voltage  prediction  for  the  25  eclipse  passages  listed  in  Table  1 (~230  potential  ob- 
servations) have  been  plotted  vs  observed  values  for  ATS-5  and  ATS-6.  The  standard 
deviation  set  is  ± 700  volts.  Error  bars  of  ± 500  V are  shown.  In  Figures  10  and  11 
four  different  eclipse  passages  are  plotted  to  give  an  idea  of  the  power  of  the  tech- 
nique. The  bars  for  the  predicted  values  correspond  to  the  eclipsed  (upper)  and 
sunlit  (lower)  predictions. 


OBSERVED  POTENTIAL  (V) 


Figure  9.  The  Observed  and  Predicted 
Potentials  for  ATS-5  and  ATS-6  for  the 
Entire  Data  Base.  The  values  for  ATS-5 
are  the  averages  of  the  sunlit  and  eclipsed 
spectra.  For  ATS-6  the  different  spectra 
are  indicated 
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Figure  10.  The  Observed  and  Predicted 
(average  of  the  sunlit  and  eclipsed  spectra) 
Potentials  for  the  Eclipse  Entry  of  ATS- 5 
on  Day  260,  1970  and  for  Eclipse  Exit  on 
Day  292,  1970 
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Figure  11.  The  Observed  and  Predicted 
(average  of  the  sunlit  and  eclipsed  spectra) 
Potentials  for  the  Eclipse  Entry  of  ATS-6 
on  Day  66,  1976  and  Eclipse  Exit  on 
Day  69.  1976 


6.  CONCLUSIONS 


The  detailed  calculation  of  the  effects  of  a time-varying  photoelectron  flux  on 
spacecraft  potential  has  hinged  on  the  concept  of  X , the  minimum  ray  path  altitude. 
Use  of  this  parameter  allowed  the  division  of  the  calculation  into  three  steps.  1-  irst, 
an  atmospheric  attenuation  function  was  derived  as  a function  of  X . Next,  follow- 
ing  the  calibration  of  the  secondary  and  backscattered  flux  ratios,  the  residual  flux 
(the  photoelectron  current)  as  a function  Xm  was  calculated  for  the  entire  set  of 
data,  f inally,  the  total  photoelectron  flux  was  determined  for  ATS-5  and  ATS-6  by 
fitting  the  atmospheric  attenuation  function  to  the  data.  Using  this  value  it  was  then 
possible  to  estimate  the  spacecraft  potential  variations  as  a function  of  time-varying 
photoelectron  flux. 

The  implications  of  the  model  are  obvious.  Tor  ATS-5  and  ATS-6  the  model 
can  predict  the  range  of  charge  buildup  as  the  spacecraft  move  in  and  out  of  eclipse 
and  provide  important  information  for  spacecraft  designers.  The  effects  of  a time- 
varying  current  source  (the  photoelectron  current)  on  spacecraft  potential  can  be 

studied.  The  model  also  evaluates  the  accuracy  of  the  simple  charge-balance 
7 9 

model  ' presented  in  Section  4.  Though  not  discussed,  the  effects  of  neglecting 
the  plasma  sheath  surrounding  the  spacecraft  can  be  evaluated.  In  the  final  analysis, 
the  major  importance  of  the  model  is  that  it  can  be  used  to  evaluate  the  interaction 
between  the  ambient  environment  and  the  spacecraft  under  commonly  occurring  con- 
ditions at  geosynchronous  orbit.  It  is  to  this  end  that  this  study  was  instigated. 


9.  Def  orest,  S.  E.  ()  97 2)  Spacecraft  charging  at  synchronous  orbit,  J.  Geophys. 

Res. . 77(No.  4):651. 
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Appendix  A 

FORTRAN  Programs 

The  FORTRAN  programs  designed  to  calculate  the  atmospheric  attenuation 
function  will  be  briefly  described.  Listings  are  provided  following  this  section.  It 
will  be  assumed  that,  as  input,  one  has  available 

RC  = angular  distance  in  radians  between  center  of  earth  and 
sun-satellite  line, 

[X  in  Eq.  (10)] 

RE  = radius  of  earth  in  R , 

G 

RS  = angular  radius  of  sun  in  radians  in  Eq.  (10)], 

XS  = distance,  in  R^,  from  satellite  to  center  of  earth 
(X S in  Eq.  (10)]  . 

The  subroutine  will  return: 

AILL  = Area  illuminated  or  percent  of  unattenuated  solar  flux 
[P  (Xm)  in  Eq.  (11)]  . 

Also  required  is  an  atmospheric  attenuation  model  as  a function  of  altitude  above 
the  earth's  surface,  X.  A model  of  the  form  given  in  Eq.  (9)  is  assumed  here.  Z0 
and  DZ  <Z(}  and  6 Z)  are  in  COMMON/ATMOS/. 


Given  RC,  RE,  RS,  and  XS,  the  subroutines  calculate  AILL  by  means  of  Eq. 
(10)  which  is  a double  integral  over  d 6 and  d <p.  The  variables  in  the  subroutine 
correspond  to  the  parameters  in  Eq.  (10)  as  follows 
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F 

4—* 

|0'  (o)  sin  1)  1 

AA 

X' 

DTH 

f-* 

do  , 

TUT 

*-4 

O . 

DR 

4 — » 

dO  . 

R 

4— ♦ 

0 . 

The  basic  program  calls  are  structured  as  follows 


N is  the  step  size.  In  the  listings  shown,  N = 10  giving  run  times  of  ~0.  02  sec/ 
value.  For  the  final  run,  N=  100  giving  run  times  of  ~ 1. 5 sec/ value. 

The  subroutine  IN  TEG,  as  indicated,  finds  1)11  and  H(I).  In  subroutine  ANGIN, 

I )T  1 1,  and  TMT(I)  are  determined.  In  subroutine  FUNG,  using  11(1)  and  TUT  (I),  F 
is  found.  Given  the  vector  TH(J),  it  is  integrated  in  QSF  to  give  THE.  THE  for 
various  values  of  R(I)  are  stored  in  HAD  (I).  HAD  (I)  is  integrated  in  QSF  to  give 
Z,  the  amount  of  the  solar  disk  visible  from  the  satellite. 

The  subroutine  QSF  is  a standard  systems  program  which  takes  a vector  Y of 
function  values  spaced  an  equal  distance  II  apart.  The  integrals  are  returned  in 
another  vector  Z. 
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SJ090UT[  f'i  rNTEG(RC,RE,RS,XS,AILLI 
OINENSU  » ?«0(100I,Z  (100) 

0R=RS/9. 

oo  l I*  ■ *1  n 

R=0R*FL  er(I-l) 

CALL  ANrT^ (R.RC.PE.XS.THE) 

1 RA0<T>=  -p- 

CALL  09  IO*, PAD, Z, 10  » 

A I LL-  7 • r ) 

»SA=f>.»  T10*  (1  .-COS(RSI  ) 

AILL^ATl '/*>A 

RETURN 

FNO 


SJOPOUT’ti-  -Ur»C(F,«C.«E,TMT,F,X$> 

conmon/  ■mos/zo.oz 

4 4 = ACOE  nC(->»*COSCRCI»SXN(RI*SIN<RCI  *C0 S(TMTt) 
X = XS*SIN«4I 
X=X-»F 
X = X*637  * . 

7 = (X-ZO  * /O* 
zisaost*' 

IC(7T.LT.1.E-G8)Z=1.E-08*Z/ZI 
IF  (71  • Of  .500 .1  Z=Z*500./ZI 
7 = EXP(-” 

z i =aos r ’> 

IF tzi.or .5  00  . ) Z=Z*500./ZI 
IF<7r.LT.l .£-08tZ-l.E-0a*Z/ZI 
7 = EXP(-7  t 
IFIX.LT.  f.  tZ-0. 

F=SIN<R)  »? 

RETURN 
FNO 

I SUOPdUT  >E  4NGIN(R,RC,RE ,XS, THEt 

niMrtlST  * TH  (1  00  ) ,7  < 1UQ  ) 

OTH=h.  7 1 ’l  0/9  . 

00  1 1=’  ,10 
T rt  T =OT  H ' ' „ 04  T ( I - 1 ) 

C4LL  F'lKM?,  PC  ,Rt  ,THT,F  ,XS) 

1 TH(Tt=P 

CALL  00  'ITN.TH.Z.IO  ) 

TME=7 (tit 
r>r  TURN 

I FNO 
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SUO’t'TTNF  osr 


(»  |pr»  • «r 

fi  iDWPUTt  THt  VICTOR  OF  INTEGRAL  VAlUT  S FOP  A GIVEN 
F '"niSTANT  TABLE  OF  FUNCTION  VALUfS. 


osar,  i 

GAI  - 1SF  <H,Y,2,NOIM) 

HE  S'*'  TP  T I 0 N OF  PARAMETERS 

H - THF  INCREMENT  OF  ARGUMENT  VAl Ut  r m 

y - THE  INPUT  VECTOR  OF  FUNCTION  VALUES. 

7 - THE  RESULTING  VECTOR  OF  INTEGRAl  VALUES.  Z MAT  BE 

I CENT  1 CAL  WITH  Y. 

N *M  - THF  DIMENSION  OF  VECTORS  Y ANO  7. 

r?C  w i\  1 ■ ~ 

N TCI  ION  IN  CASE  NOIM  LESS  THAN  3. 

SUB’’ H TTMES  ANO  FUNCTION  StJOPPOGFAMS  kEQUIFFO 

•IO»r 


Mr  TM  " 

T '’INNING  WITH  ZI1>=0,  IVAUJATION  OF  VECTOR  7 TS  OONE  BY 
M 'V 5 Of  SIMPS  INS  RULE  TOGETHER  WITH  NFWTON'’  3/P  RULE  OP  A 
~ " TTNATION  OF  THT  St  TWO  RULES.  TRUNCATION  reerjp  is  0F 

0 H**S  II . E • FOURTH  ORDER  METHOD).  ONE v TN  CASE  NOIM=J 
T INCATION  ERROR  OF  ZI2)  IS  DF  ORDER  H**4. 

F f ?F  Ff  RE  NOE,  SEE 

1 I F.t.HILUEURANO,  INTPOOUCTION  TO  NUMERICAL  ANALYSIS. 

h:gfaw-hill,  new  yofk/toronto/london,  idgf,  pp.n-zf,. 

I I P.  ZURHIJLHL  , PRAKTISCHr  M A THE  HA  T I K FUFP  IMF,  r N IFUPE  UNO 
RR  Y S I K f P . SPRINGER,  DF  RL  I N/ GCE  T T INGE  N/ HFIOFL  BERG,  1963, 
DJ .214-221 . 


SUTPOUrt'  irf  <H,Y,Z,NOIM) 


OIMrrrT''*  Y(1  I ,z  ID 
HTs.333’  w 3*H 

i f i >m  m - r 1 7 , h , i 

N’)IM  IS  C?-ATIP  THAN  b.  PREPAPATIDNS  OF  INTEGRACTDN  LOOP 
t SUM  1 = Y I '*YI?| 

SUN  1 - S U USJM1 

SJH1=HT  • IM  n ♦ SUM1TY  I 3»  I 

A JXlsYI't  *»VI4I 

AiKlsAU  *!♦  A J XI 

AUX1  = SUMDHT*  IVI3)*AUX1*YI5)  I 

A'JX?=HT»  IYII  ) ♦3.BP5*  I Yl  21  AVI  DM  IYI  3)*Y  I4M*V|BI  » 

SUN?=Y  I 5 *♦  Y I *!  | 

'•l)N?=SU1  "♦  SJM2 

S )M?=AU  ’-NT* I Y I 4 ) ♦ $ UM2 ♦ Y I b ) ) 


-A. 
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7(11=0. 

4IJX=Y(3>  *f  (3) 

AUX=AUX*  MX 

Z(2)=SUN?-HT* (Y(21*AUX*Y(4) ( 

7 ( 31  =SU  I 
7( 4>=SUt  ? 
rF(W)IM-Y»0,5,2 

INTFC.R4' TON  LOOP 

2 00  4 1*7  ,NOI  M,  2 
S'J*i  1= a i ’i 
SUM?=  All'  ? 

A'lX  1 = V »'  -1  >*Y  (1-1) 

AUX1  = AljY  U 4UX1 

A'JXl*SU  M»HT*  ( Y<  1-2)  ♦AUXl  + YJ  III 
7 ( I -2)  all'll 
TFII-NTT' I 3,8,6 

3 A'JX7  = Y ( ’»  Y(  I» 

AUX2=A'J  v?»  4 JX2 

4UX?=SU*  >HT*  (Y(  1-1 1 YAU X2*Y(  I*1M 

4 7(1-11  «V.1 1 2 
1 7 ( NOI M-1*=4JX1 

7 ( NOTMI  = " JX2 
RR TURN 

6 7 ( NOI M- 1 TrOJM? 

7 ( NOT  M I rAjXl 
Rr TURN 

f NO  OF  T*K  Fi  J/i  T I ON  LOOP 

7 IF  INOU-  -0  12 , 1 1,  6 
N 1 1 M 10  n'»L  TO  4 OR  5 

8 TIM ?=  1 . 1 r5»HT»(Y(l>*Y(2)»Y(2IYY(2)»Y(3>»Y(3l»Y(31*Y(41l 
F'JNlrV  ( ’OTI2I 
3'JMlaSU  1 •*  0JH1 
S'H1  = MT*  If  (1  ) ♦SUM1»Y  ( 3»  I 
7 I 1 1 = 0 . 

AUX  1=Y  f 3 >♦  Y ( 3 ) 

A JX1=4UV1* 4JX1 

7(,)  = 3U'  2-HT  • ( Y(  2>»AUX1*Y(4>  I 
IF(N0I1-rl  10,9,1 

9 4 'J  X 1 = V ( 1 ’»  Y < 4 I 
A'JXt=4'l7  1 » 1UX1 

7(r.l  = 9U’  U (I*(Y(3»*flUXl»Y(5)  I 

10  7 ( 31 =3U ’ * 

7(4l  = Vr  ' 

O - TUDN 

C 

C N T I M 10  TTJ4L  TO  3 

11  0 1 J N 1 n M T M .2*»Y(  1)»Y  (2)  ♦Y(2I -,?9*Y(3)  ) 

O'M 2-YC  *fY(?l 

S'JM?=5U  ' "f  3J"2 

7 ( 3 I = H T II  (1  ) ♦ SU  M 2 ♦ Y ( 31  I 

7(11-1). 

7 ( ->|i5U  1 1 

12  RFTU7N 
FNO 
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